Improved Time-Domain Accuracy Standards for Model Gravitational Waveforms 
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Model gravitational waveforms must be accurate enough to be useful for detection of signals and 
measurement of their parameters, so appropriate accuracy standards are needed. Yet these standards 
should not be unnecessarily restrictive, making them impractical for the numerical and analytical 
modelers to meet. The work of Lindblom, Owen, and Brown [Phys. Rev. D 78, 124020 (2008)] 
is extended by deriving new waveform accuracy standards which are significantly less restrictive 
while still ensuring the quality needed for gravitational-wave data analysis. These new standards 
are formulated as bounds on certain norms of the time-domain waveform errors, which makes it 
possible to enforce them in situations where frequency-domain errors may be difficult or impossible 
to estimate reliably. These standards are less restrictive by about a factor of 20 than the previously 
published time-domain standards for detection, and up to a factor of 60 for measurement. These 
new standards should therefore be much easier to use effectively. 
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I. INTRODUCTION 



Expected astrophysical sources of observable gravtia- 
tional waves are often systems dominated by strong grav- 
itational dynamics. In many cases this allows clear de- 
tailed predictions for the gravitational waveform signals, 
based on Einstein's theory, General Relativity While 
these theoretical predictions are well-defined in principle, 
explicit computation of the model waveforms requires a 
combination of challenging analytic [l[ and numerical (H 
computations to approximate the ideal theoretical expec- 
tations. If the approximation errors are too large they 
will impact the interpretation of gravitational wave data 
from operating gravitational wave instruments such as 
LIGOjp and Virgo 0, or future instruments such as 
LISA [| and LCGT [J. 

Model waveforms are used for matched filtering in 
gravitational-wave data analysis both for detection and 
for measurement of the physical parameters of the source 
(referred to as parameter estimation by the data analy- 
sis community). A signal is first identified in a detector's 
noisy data stream when it is found to have a significantly 
large correlation with (i.e. a sufficiently large projection 
onto) a model waveform. If the model waveforms were 
not accurate enough, then an unacceptably large fraction 
of real signals would fail to be detected in this way. The 
second use of model waveforms is to measure the phys- 
ical properties of any signals identified in the detection 
step. These measurements are performed by fine tun- 
ing the model- waveform parameters (e.g., the masses and 
spins of the components of a compact binary, the source's 
orientation, the times of arrival of the signals, etc.) to 
achieve the largest correlation with the data. If the model 
waveforms were not accurate enough, these measured pa- 
rameters would fail to represent the true physical proper- 
ties of the sources to the level of precision commensurate 



with the intrinsic quality of the data. The systematic 
errors due to modeling inaccuracies should be no greater 
than the statistical errors due to detector noise. Model 
waveforms produced by numerical relativity simulations, 
e.g. of binary black holes, are computationally intensive 
to generate, and thus it is desirable to know what the 
minimum accuracy standards needed for these two uses 
are. Flanagan and Hughes Q made the first derivation 
of abstract and generic accuracy standards in this con- 
text. Miller Q made the first concrete application of 
these standards to binary black-hole waveforms in nu- 
merical relativity, comparing post-Newtonian errors to 
then-current simulations and arguing for improvements 
in the latter. Our recent work has focused on finding 
simple intuitive derivations of the fundamental accuracy 
standards, extending the standards to include the effects 
of detector calibration errors, and deriving new represen- 
tations of the standards that are easier for the numeri- 
cal relativity community to verify the waveform accuracy 
needed for gravitational- wave data analysis (9j-[Tl| . 

These accuracy standards are expressed as limits on 
the waveform- modeling errors Sh m = h m — h e , the differ- 
ence between a model waveform, h m , and its exact coun- 
terpart, h e . Numerically modeled gravitational wave- 
forms, and the errors associated with them, are gener- 
ally evaluated most easily as functions of time, h m (t) 
and Sh m (t). In contrast gravitational- wave data analy- 
sis and the accuracy standards for model waveforms are 
most conveniently formulated in terms of the frequency- 
domain representations of the waveforms and their errors, 
hm(f) and Sh m (f). The time- and frequency-domain 
representations are related to one another by Fourier 
transforms, e.g., 
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h m (t)e- 2mft dt. 



(1) 



So it is straightforward (in principle) to transform from 
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one representation to the other. 

The simplest way to express the standards needed to 
ensure the appropriate levels of accuracy for model grav- 
itational waveforms is to write them in terms of a par- 
ticular norm of the model- waveform errors: (5h m \5h m ) . 
This norm, defined by 



(Sh m \Sh r , 



6h m (f)Sh* m (f) 

Sn(f) 



df, 



(2) 



weights the different frequency components of the wave- 
form error by the power spectral density of the detector 
noise S„ (f). In terms of this norm, the accuracy require- 
ment that ensures no loss of scientific information during 
the measurement process is, 



(Sh m \Sh m ) rjc 
(h m \h m ) p 



(3) 



cf. See s. VI and VII of Ref. and Eq. (5) of Ref. 0. Here 
P = v (h e \h e ) is the optimal signal-to-noise ratio of the 
detected signal, and the parameter r\ c < 1 is set by the 
level of calibration error in the detector, cf. Appendix |A1 
and Ref. [lOj . Similarly the accuracy requirement that 
ensures no significant reduction in the rate of detections 
is, 



(Sh m \Sh m 
(h m \h m ) 



< V2e 



(4) 



where e m ax is a parameter which determines the fraction 
of detections that are lost due to waveform-modeling er- 
rors, cf. Eq. (14) of Ref. @. The choice e max = 0.005 en- 
sures that no more than 10% of the signals will be missed, 
assuming the template bank discretization currently used 
in LIGO searches for coalescing compact binaries [9j. 

The basic accuracy standards, Eqs. Q and give 
the needed bounds on the waveform errors Sh m . Un- 
fortunately, these abstract requirements are difficult, or 
perhaps impossible, to enforce directly in practice [111 ]. 
One fundamental problem is that an exact knowledge of 
the error, 8h m , associated with a model waveform, h m , is 
never known. If Sh rn were known, then the exact wave- 
form h e — h m ~ Sh m would also be known, and there 
would be no need for accuracy standards. At best, the 
waveform modeling community can aspire to construct 
tight upper bounds on these errors, \Sh m \ < SH m , that 
would make it possible to guarantee the needed waveform 
accuracy standards. Such bounds are often very difficult 
to construct. So in the absence of such bounds, good 
estimates of the errors, 6h m ss SH rn , are a good way to 
verify that the accuracy standards are satisfied at least 
approximately. 

Another (more practical) problem with the basic accu- 
racy standards, is that they are formulated in terms of the 
frequency-domain errors Sh m (f). Since numerical model 
waveforms are generally computed in the time domain, 
h m (t), it will almost always be easier and more straight- 
forward to bound or estimate the time-domain waveform 



errors 5h m (t). Time-domain quantities are converted to 
the frequency-domain via the Fourier transform, Eq. ([1} . 
In principle this is straightforward. In practice, how- 
ever, it is easy to introduce errors in this process that 
can be larger than the intrinsic waveform modeling er- 
rors. Discontinuities and other non-smoothness at the 
beginning and end of a finite length waveform can cre- 
ate significant errors in the frequency-domain represen- 
tation, unless appropriate "windowing"_procedures are 
followed. This can be done (cf. Ref. [12| for a recent 
suggestion) but it requires some care and has some lim- 
itations on applicability. Also, more fundamentally, a 
time-domain bound does not translate to a frequency- 
domain bound [ll|: that is, Sh m (t) < SH rn (t) does not 
guarantee Sh m (f) < SH m (f). For reasons of convenience 
and reliability therefore, it is very desirable (perhaps even 
necessary in some cases) to have versions of the accu- 
racy standards formulated directly in terms of the time- 
domain waveform errors Sh rn (t) . So this paper is devoted 
to finding time-domain representations of these standards 
that are straightforward and practical to use without be- 
ing overly restrictive. 

The most straightforward approach to obtaining the 
needed time-domain expression for the waveform accu- 
racy would be to re- write the fundamental noise weighted 
norm (5h m \5h m ) in terms of a suitably defined time- 
domain inner product. This can be done exactly us- 
ing the "noise kernel" k(i), defined as the real part of 
the Fourier transform of 2S~ 1 (f). The resulting time- 
domain norm is then identical to the standard noise 
weighted frequency-domain norm (e.g. Ref. (l3j): 

oo poo 



(Sh rn \Sh m ) = / 5h m (ti)k(ti - t 2 )5h m (t2)dtidt 2 . 

J —oo J —oo 

(5) 

Since this norm is identical to the standard noise- 
weighted norm, it can be substituted directly into the 
fundamental accuracy standards, Eqs. © and @. Un- 
fortunately for the same reason, it also suffers from both 
of the problems discussed above. Upper bounds in the 
time domain |5/i m (i)| < 6H m (t) do not produce upper 
bounds on these time-domain norms. And errors in the 
norms introduced by Fourier transforms of finite dura- 
tion waveforms with non-smooth edges (Gibbs phenom- 
ena, etc.) also occur, because this time-domain norm is 
non-local in time on the timescales to which the detector 
is sensitive. These problems could be avoided if there 
were efficient representations of the accuracy standards 
that placed limits on local-in-time norms of the waveform 
errors. The purpose of this paper is to construct new for- 
mulations of the waveform accuracy standards that meet 
these criteria. 

Before moving on to a detailed discussion of time- 
domain error standards themselves, it is worth pausing 
to think briefly about how waveform errors 5h m have 
been (or could be) bounded or estimated. Error bounds 
(called a posteriori error estimates in the mathemat- 
ics literature) have been constructed for quantities de- 
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rived from the numerical solutions of hyperbolic evolu- 
tion problems [l4j]. However no results of this kind are 
known at this time (to our knowledge) specifically for 
the case of gravitational waveforms extracted from so- 
lutions of Einstein's equations. The error bounds that 
do exist are for time-domain errors, and we know of no 
frequency-domain bounds for quantities derived from the 
solutions of time-evolution problems at all. Given a time- 
domain error bound, SH m (t) , it would only be useful for 
enforcing the basic accuracy standards if it could be con- 
verted to a frequency domain bound \Sh m (f)\ < 5H m (f), 
or if a time-domain version of the accuracy standards 
were available. Unfortunately the Fourier transform of 
5H m (t) does not generally provide the needed frequency- 
domain bound, since (SH m \8H m ) is not always larger 
than (Sh m \Sh m ) (llj . Time-domain bounds, 5H m (i), are 
only useful therefore if there are versions of the basic ac- 
curacy standards, Eqs. (J3|) and (gj, based on local-in-time 
norms of the waveform error. 

Good approximations of the waveform errors Sh m (t) ~ 
6H m , can be constructed using Richardson extrapola- 
tion methods. Model waveforms generally have errors 
that scale with some accuracy parameter, a, in a well 
understood way. 1 Many methods of computing wave- 
forms have dominant errors that scale as some power of 
this parameter, Sh m (t) ~ A(t)a n , while others have er- 
rors that scale exponentially, Sh m (t) ~ A(t)e~ a ^ a . In ei- 
ther case the values of the quantities A(t), n, a, etc. can 
be evaluated by comparing the waveforms h m (t, aj) pro- 
duced by simulations using different values of a,;. The er- 
rors in the most accurate of these simulations, a = a m - m , 
is then given approximately by the scaling expression, 
e.g., Sh m (t) « SH m (t) = A(t)a™ in , using the measured 
values of the various constants. These Richardson ex- 
trapolation methods could also be applied to frequency- 
domain representations of the waveforms evaluated at 
different accuracies, h m (f,a,i) (subject to the problems 
involved in performing the Fourier transforms discussed 
above). The resulting frequency domain error estimates, 
e.g., Sh m (f) « SH m (f) = A(/)a™ in , could then be used 
directly to estimate whether the fundamental accuracy 
standards are satisfied. The existence of accuracy stan- 
dards based on local-in-time norms just provides a more 
convenient and perhaps a more reliable way to evalu- 
ate waveform accuracy in this case, rather than being a 
necessity as it was in the error-bound case. The gravita- 
tional waveform simulation community has begun to es- 
timate waveform errors by comparing time-domain wave- 
forms computed with different accuracy parameters [l5r - 
[22j . However, Richardson extrapolation methods have 
not yet been applied to improve the quality of these er- 
ror estimates (to our knowledge). 



The discussion above shows that alternate versions 
of the basic accuracy standards, based on local-in-time 
norms of the waveform errors, arc needed for reasons 
of convenience and perhaps necessity. This paper con- 
structs a number of new versions of the standards that 
meet all the needed criteria. Some representations of the 
accuracy standards based on local-in-time norms of the 
time-domain errors have already been discussed in the 
literature [t| [ll| . While these are sufficient to guaran- 
tee the fundamental accuracy requirements, Eqs. §3§ and 
(U]) , they achieve this at the price of placing excessive re- 
strictions on the allowed time-domain errors. New time- 
domain accuracy standards are developed here that are 
considerably less restrictive: about a factor of 20 for the 
detection standards, and up to a factor of 60 for mea- 
surement. These new standards are at most a factor 
of three more restrictive than their optimal frequency- 
domain counterparts, so we expect they should be prac- 
tical to use. The previous work on time-domain accuracy 
standards is reviewed, and the new standards are derived 
and discussed in detail in the following sections of this 
paper. 



II. TIME-DOMAIN ACCURACY STANDARDS 

The fundamental frequency-domain accuracy stan- 
dards, Eqs. © and (@|, can also be expressed in terms of 
the time-domain L 2 norm of the waveform errors 0, [ll[ . 
These time-domain versions provide a way to enforce the 
waveform-accuracy standards without the need for a de- 
tailed knowledge of the frequency-domain representation 
of the errors. This section reviews the derivation of the 
presently known time-domain standards, and illustrates 
the shortcomings that make them an impractical way 
to enforce the accuracy requirements. The methods of 
analysis used to derive the original time-domain stan- 
dards arc then generalized to produce several new and 
improved versions. Each of these new standards is itself a 
sufficient condition that can be used to enforce the funda- 
mental frequency-domain standards. As an illustration, 
these new standards are applied to the case of binary 
black-hole inspiral-merger-ringdown waveforms for a de- 
tector with an Advanced LIGO noise curve 



23J. This 



example shows the effectiveness and utility of these new 
time-domain standards as tools for waveform-accuracy 
enforcement. 

The frequency-domain accuracy standards, which are 
based on the noise- weighted norm (Sh rn \Sh m ), can be 
converted to time-domain standards using the following 
basic inequality: 



(Sh m \6h r . 



< 



4f™\6h m (f)\ 2 df 
min[5„(/)] 



(6) 



This parameter could represent the spacing between grid 
points for numerical waveform calculations, or perhaps a post- 
Newtonian like expansion parameter for an analytic calculation. 



This inequality approximates the power spectral den- 
sity of the detector noise as a constant: its minimum 
value. The approximation in Eq. © is very good there- 
fore whenever the model waveform error is largest in the 
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sensitive band of the detector, and is not so good when 
a significant part of the error lies outside this sensitive 
band. By Parseval's theorem, the numerator in Eq. ([6]) 
is proportional to the time-domain L 2 norm, ||<5/i m (i)||: 



ll*M*)ll 



2 _ 



\5h m {t)\ 2 dt = 2 / \6h m {f)\ z df. (7) 

c JO 



Thus the basic inequality, Eq. ©, can be re- written as 



(Sh m \Sh r 



< 



2\\Sh m (t)\\ 2 
min[5„(/)] 



(8) 



Using this inequality, the fundamental accuracy stan- 
dards, Eqs. and ((H), can be re- written in terms of 
time-domain L 2 norms: 



\\Sh m (t)\\ <Co r k^ 
\\h m {t)\\ ~ p 1 



for measurement, and 

\\Sh m (t)\\ 

IIM*)II 



< CoVTe 



(9) 



(10) 



for detection. The quantity Co that appears in these ex- 
pressions is defined as the ratio of the traditional optimal 
signal-to-noise measure p to a non-traditional measure: 



r 2 - n 2 



2\\hm(t)\\ 2 
min[5„(/)] 



(11) 



This quantity is dimcnsionless and independent of the 
overall scale (i.e. the distance to the source) used in h m . 

The time-domain standards, Eqs. © and (fTU)) , differ 
therefore from the fundamental standards, Eqs. ([3]) and 
in just two ways: i) the norms used to measure the 
waveform errors on the left sides are different, and ii) 
the maximum allowed errors on the right sides are mul- 
tiplied by the factor Cq in the time-domain case. It is 
straightforward to show that Co is always less than one: 
Cq < 1. Thus the time-domain standards are always 
more restrictive than their frequency-domain counter- 
parts. The quantity Co compensates in the time-domain 
standards for the fact that the waveform errors are not 
being weighted by the detector noise in the optimal way. 
Unfortunately these time-domain standards have two se- 
rious flaws: The first is that the time-domain norms 
||/i e (i)|| are not well defined, because the waveforms h e (t) 
do not fall to zero quickly enough as t — > —oo. Since most 
model waveforms h m (t) are only evaluated for finite time 
intervals, this problem is one of principle more than prac- 
tice. But this does mean that the norms ||/i m (t)||, and 
thus the quantity Co, depend (although fairly weakly) on 
the length in time of the model waveform. The second 
flaw, discussed in more detail below, is that Co is quite 
small for realistic waveforms, thus making these time- 
domain versions of the standards very over-restrictive. 

To determine how overly restrictive the time-domain 
standards are in practice, the quantity Co has been evalu- 
ated for non-spinning equal-mass black-hole binary wave- 
forms and for a detector with an Advanced LIGO noise 
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FIG. 1: Curve illustrates Co, the ratio of the traditional op- 
timal signal-to-noise measure p to a non-traditional measure 
defined in Eq. {TTJ), as a function of the total mass for non- 
spinning equal-mass black-hole binary waveforms and for a 
detector with an Advanced LIGO noise curve optimized for 
double neutron star binaries. 



curve optimized for double neutron star binaries [231] . 2 
These waveforms were constructed by patching together 
post-Newtonian waveforms for the early inspiral phase, 
with numerical relativity waveforms for the late inspiral 
plus merger plus ringdown phases (2~H . |25| and adjust- 
ing the time scales for different total masses. They are 
the same used earlier to evaluate Co . Figure Q] il- 
lustrates Co as a function of mass for binary black-hole 
systems with masses in the range that is relevant for Ad- 
vanced LIGO: 4 < M/Mq < 400. For low mass bi- 
nary systems — where these model waveforms lie mostly 
within the Advanced LIGO frequency band — Co is of or- 
der unity, and the time-domain standards arc reasonably 
efficient representations of the fundamental frequency- 
domain standards. However, for binary systems with 
large masses — where these model waveforms lie mostly 
outside the Advanced LIGO band — Co becomes very 
small, and the time-domain standards become very inef- 
fective. For the largest masses the quantity Co becomes 
very small, Co ~ 0.016, and the time-domain standards 
must place unreasonably tight constraints on the wave- 
form errors to ensure that the fundamental frequency- 
domain standards are satisfied by the small part of these 
waveforms lying within the Advanced LIGO band. 



2 The data for this neutron-star optimized noise curve comes from 
the file NSNS.Opt.txt available from the URL in Ref. [H- 
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The methods of analysis used to derive Eqs. (|9]) and 
(fTUj) can be generalized to obtain a new set of time- 
domain accuracy standards. This can be done using the 
following extension of the basic inequality: 



(8h m \8h Tl 



< 



Aj™(2nf) 2k \5h m {f)\ 2 df 
min[(27r/)2^„(/)] ' 



(12) 



The numerator in this expression is just the frequency- 
domain L 2 norm of the fc th time derivative of the wave- 
form error. This can be converted to a time-domain L 2 
norm using Parseval's theorem: 



l<*M*)llfc 



d k Sh m (t) 



dt. 



= 2 



dt k 

(2irf) 2k \5h m (f)\ 2 df. (13) 



The basic inequality in Eq. (jl2j) can therefore be re- 
written as 



(5h m \6h r , 



< 



min[(27r f) 2k S n (f)Y 



(14) 



These inequalities together with the fundamental 
frequency-domain accuracy standards, Eqs. ([3]) and Q, 
can then be used to obtain a new set of time-domain 
accuracy standards: 



£k < Cfc — , 

p 

for measurement, and 

£k < Cfc\/2e max , 
for detection. The quantities £k, 

k ~ \\h m {t)\\k ' 



(15) 



(16) 



(17) 



are new measures of the time-domain waveform error, 
and the quantities Ck are defined as the ratios of the 
traditional optimal signal-to-noisc measure p to new non- 
traditional measures: 



C 2 - i? 



2|lM*)lli 



min[(27r/) 2 ^ n (/) 



(18) 



Both quantities, £k and C^, are dimensionlcss and in- 
dependent of the overall scale (i.e. the distance to the 
source) used in h m . It is straightforward to show that 
Cfc < 1, so these new standards are always more re- 
strictive than their frequency-domain counterparts. The 
k = versions of Eqs. (fT5|) and (fT6|) arc identical to the 
previously known time-domain standards, Eqs. © and 
(fTOj) . But the standards with k > 1 are new and distinct. 
The new standards with k > 1 arc free from the flaws 
described earlier for the previously known fc = stan- 
dards. The norms ||/i m (t)||fc are well defined for fc > 1, 
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FIG. 2: Curves illustrate Ck, the ratio of the traditional op- 
timal signal-to-noise measure p to a non-traditional measure 
defined in Eq. (|18|) . as a function of the total mass for non- 
spinning equal-mass black-hole binary waveforms. The wave- 
form accuracy standards, Eqs. (|15[) - (|16[) , can be enforced us- 
ing any value of k. Thus it is sufficient to impose the standard 
having the largest Ck in each mass range. 



because d k h m /dt k falls off quickly enough as t —> — oo. 
And as described in more detail below, the Ck for fc > 1 
are much closer to one. So these new standards are much 
more effective ways to enforce the accuracy requirements. 

Figure [2] illustrates the quantities Ck (for fc = 0,1, 
and 2) that appear in the waveform accuracy standards, 
Eqs. (fT"5j) and (fT6|) . as a function of the mass of the wave- 
form. Each of these new accuracy standards is sufficient 
to guarantee that the corresponding frequency-domain 
standard is satisfied. Choosing the easiest standard to 
satisfy is all that is required for a given M, and this will 
typically be the one having the largest Ck ■ Figure[2]shows 
that the standards with different values of fc have differ- 
ent mass ranges where they are most effective. The fc = 
standard is most effective for low-mass binary waveforms, 
while the fc = 2 standard is the most effective for high- 
mass waveforms. 

Figure [3] shows the noise curves used to compute the 
C k illustrated in Fig. EI These particular noise curves 
are all representations of the Advanced LIGO sensitiv- 
ity curves optimized for neutron-star/neutron-star binary 
signals [23j. The noise curves with successively larger 
values of fc have minima at successively smaller values of 
the frequency /. Figure [3] shows that the factor f 2k com- 
bines with the piecewise power law structure of S n (f) 
to produce fairly broad regions where the effective noise 
curves are rather flat. In these flat regions the approx- 
imation, min[/ 2fc 5' n (/)] ~ f 2k S n {f), becomes extremely 



6 



10 



10 



10 



10 



10 



-24 



30 



36 



42 



48 



1 1 


i 




(2 7tf) S n (f) 




k = 




- 


— - k = 1 






k = 2 




4-. 

i •••• 
1 
















•— 1 ^^^^ 


■ i 


i 


1 



io 1 



10 10 
f (Hz) 



io H 



FIG. 3: Curves illustrate the effective noise curves, 
(2irf) 2k S n (f), that appear in the norms of the k time 
derivatives of the waveform. The minima of these curves are 
marked by the large dots. These curves illustrate why the 
higher time derivative norms are more sensitive to the lower 
frequency parts of the waveforms, and hence do a better job 
of modeling the waveform accuracy in high mass binaries. 



good, and the approximation leading to Eq. (fI2)l be- 
comes nearly exact. These features explain one of the 
reasons the higher time-derivative norms measure the 
lower-frequency parts of the waveforms more faithfully, 
and hence do a better job of enforcing the waveform ac- 
curacy standards in high-mass waveforms. The higher 
time-derivative norms measure d k h/dt k , which is more 
strongly peaked near the time of merger. Hence the 
energy norms Ek for k > 1 do a better job than Eq of 
measuring the strongest parts of the waveform that con- 
tribute most to detection. 

Model waveforms used to measure accurately the phys- 
ical parameters of a previously detected signal may be 
custom made for this task. The physical parameters 
of the waveform will be known (approximately) in this 
case from the detection process, so a finely targeted tem- 
plate bank of waveforms with parameters that cover a 
small neighborhood surrounding the actual physical val- 
ues can be computed. For such applications, the quan- 
tities Cfc can easily be computed along with the model 
waveforms. The waveform accuracy standard for mea- 
surement, Eq. (|15p . should be enforced therefore using 
the standard having the largest Ck for these particular 
waveforms. As Fig. [5] shows, this will typically be the 
k = 2 standard for the largest mass binary waveforms, 
k = for the smallest mass waveforms, and k = 1 for a 
broad range of intermediate mass waveforms. The small- 
est value of Ck needed for these measurement tests is 
Co = C\ w 0.331 if the mass of the waveform is near 
IOMq. However, if the mass is larger, the value of the 



Ck could be as large as C2 ~ 0.935 for waveforms with 
masses near 400 Af©. These new time-domain accuracy 
standards for measurement can therefore be up to 60 
times less restrictive (depending on the mass of the wave- 
form) than the previously available k = standard Q . 

The situation is somewhat more complicated for wave- 
forms used as part of the detection process. These wave- 
forms should be scalable to the full range of masses whose 
waveforms lie within the sensitivity band of the detector. 
The simplest way to do this would be to use the ac- 
curacy standard having the largest miiiM Ck(M), where 
the minimization is done over the waveform masses in 
the appropriate range. Note that the quantity Ek , which 
appears on the left side of Eq. ([16]) , is independent of 
mass because both ||£/i m ||fc and ||/i m ||fc scale with mass 
in exactly the same way. Thus the accuracy standard 
for a given k can be enforced over the entire range of 
masses by ensuring that Ek satisfies the standard with 
Ck = nrinM C/-(Af). The k = 1 standard has the largest 
minM Cfe(M) for the waveforms illustrated in Fig.[2J so it 
is sufficient in this case to enforce the k = 1 standard with 
C\ = minjvf C\(M) as 0.239. This is about 15 times less 
restrictive than would be required using the previously 
available k = standard where minM Co(M) « 0.016. 

A less restrictive accuracy standard for detection can 
also be constructed, at the expense of making it slightly 
more complicated. Figure [2] shows that enforcing the 
detection standard, Eq. (frp|) . with C\ = minju C\(M) is 
somewhat more restrictive than is really necessary. All 
that is required is to enforce the standard having the 
largest Cfc(M) for each mass. So it would be sufficient to 
enforce the k = standard in the mass range where Co 
is the largest of the Ck'- M < M- lnt . The mass M- lnt ss 
10.4M© represents the point where the Cq(M) and the 
Ci(M) curves intersect: C (Af int ) = Ci(A'/ int ) « 0.331 
for the waveform example shown in Fig. [2] It would also 
be sufficient to enforce the k = 1 standard in the mass 
range M > M- lnt . Therefore it is sufficient over the full 
range of masses to enforce both the k = and the k = 1 
versions of the detection accuracy standard, Eq. (|16|) . 
with C (A/ int ) = Ci(Af int ) » 0.331: 3 



E a < 0.331 y/2T„ 



and Ex < 0.331 



(19) 



This condition is about 20 times less restrictive than 
would be required to enforce the previously available 
k = standard over the full range of masses. The overall 
factor that appears on the right sides in Eq. (fTT)| appears 
to be quite insensitive to the detector noise curve, giving 
0.351 instead of 0.331 using the "zero detuning, high- 
power" version of the Advanced LIGO noise curve. 4 We 



3 In principle it would be sufficient to enforce the k = 2 instead of 
the k = 1 standard in the mass range M > 26OM0. However, 
this would not make the needed k = 1 condition any weaker, 
because S\ is independent of mass. 

4 The data for this zero-detuning, high power noise curve comes 
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have not yet determined how sensitive this factor is to 
properties of the waveform models, like the total dura- 
tion of the waveform, or the spins of the individual black 
holes. 



III. DISCUSSION 

New time-domain representations of the waveform ac- 
curacy standards for gravitational-wave data analysis 
have been presented in Sec. [TTJ For waveforms to be 
used for measuring the physical properties of a detected 
signal, the new time-domain standards, Eq. (|15j) . pro- 
vide the best way to enforce the needed accuracy stan- 
dards. These standards give the allowed upper limits 
on certain new measures of the waveform error de- 
fined in terms of the time-domain L 2 norms of the fc th 
time derivatives of the model waveform h m and its er- 
ror 5h m : E k = \\5h m (t)\\ k /\\h m (t)\\ k . Model waveforms 
are generally obtained by performing time integrals of 
an equation like d 2 h m (t)/dt 2 = ^(i), where ^(i) is the 
(real) projection of the outgoing component of the Weyl 
curvature representing the wave polarization of interest. 
Therefore, it is no harder or less accurate to evaluate the 
first two time derivatives of h m (t), and their accuracy, 
than it is to evaluate h m (t) and its accuracy. For the 
binary black-hole waveforms used as an example here; 
the k = version of this standard is most effective for 
black- hole systems with masses below about 10A/©; the 
k = 1 standard is best in the intermediate mass range 
10 < M/Mq < 260; and the k = 2 standard is best 
for the highest mass waveforms, 260 < M/M@ < 400. 
This new version of the accuracy standard for measure- 
ment is less restrictive than the previously available time- 
domain standard by up to a factor of 60 for the largest 
mass binary systems, and is never more restrictive than 
the ideal frequency-domain standards by a factor greater 
than about 3. The new standard for measurement is 
nevertheless sufficient to ensure that the ideal frequency- 
domain standard is satisfied as well. 

The situation is slightly more complicated for the new 
time-domain waveform accuracy standards for detection. 
These standards put limits on both the k = and 
k = 1 versions of the time-domain error measures, £q 
and E\. The best new time-domain standard for detec- 
tion presented in Sec. |TT] enforces these standards with 
Co = C\ w 0.331 for the example waveform used here. 
This new standard for detection, Eq. (fT9]) , guarantees 
that the frequency-domain standards are enforced, and it 
is only about a factor of three more restrictive than those 
ideal standards. The new standards should be effective 
enough therefore to make them usable without placing 
an undue burden on the waveform-simulation commu- 



from the file ZERO_DET_high_P.txt available from the URL in 
Ref. [H. 



nity. Slightly less restrictive time-domain standards for 
detections can also be derived using methods similar to 
those used in Sec. Ql] The derivation of these additional 
alternative standards is discussed and analyzed in detail 
in Appendix [B] These alternative standards are more 
complicated, however, and consequently may not prove 
to be as useful as Eq. (Til)]) , 

Acknowledgments 

We thank Michael Hoist for valuable discussions about 
rigorous mathematical error bounds, Sean McWilliams 
for other useful discussions, and Duncan Brown for help- 
ful comments on an earlier draft of this paper. This re- 
search was supported in part by grants to Caltech from 
the Sherman Fairchild Foundation, NSF grants DMS- 
0553302, PHY-0601459, and PHY-0652995, and NASA 
grant NNX09AF97G; by NASA grants 08-ATFP08-0126 
and 09-ATP09-0136 to Goddard Space Flight Center; by 
NSF grant PHY-0855589 to Penn State; and by the LIGO 
Visitors Program. LIGO was constructed by the Califor- 
nia Institute of Technology and Massachusetts Institute 
of Technology with funding from the National Science 
Foundation and operates under cooperative agreement 
PHY-0757058 . This paper has LIGO Document Number 
|ligo-pl000078 -v2. 

Appendix A: Calibration Error Levels 

The basic accuracy requirement for measurement, 
Eq. reduces to the case of an ideal detector when 
?/ c = 1. The ideal detector approximation means that 
the systematic errors made in calibrating the detector 
(by measuring the response function of the detector) are 
negligible compared to the statistical errors due to the 
intrinsic noise level of the detector. It is more realistic 
to expect that in many cases the detector will have non- 
negligible calibration errors. These errors influence the 
data analysis process in almost the same way as waveform 
modeling errors [Ioj |. Scientific information will be lost 
unless the combined systematic calibration and waveform 
modeling errors are smaller than the intrinsic detector 
noise level. This requires that both the calibration error 
and the waveform modeling errors be kept below the in- 
trinsic noise level of the detector. The allowed error bud- 
get may however be apportioned arbitrarily between the 
two sources, e.g., in the most economical or convenient 
way. The parameter r\ c introduced in Eq. determines 
the share of this error budget allowed for waveform mod- 
cling error. In particular it represents the amount by 
which the waveform modeling error must be reduced be- 
low the intrinsic detector noise level of the detector to 
allow for the presence of calibration error. 

In some cases detectors may have calibration errors 
that arc negligible compared to the needed waveform 
modeling errors (e.g., the proposed Laser Interferometer 
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Space Antenna LISA). In those cases the ideal-detector 
conditions apply and r/ c = 1 is the appropriate choice. 
In most cases, however, it may only be possible or eco- 
nomically feasible to calibrate to an accuracy level com- 
parable to the intrinsic detector noise. In these cases 
it seems likely that it will be optimal to split the avail- 
able error budget equally between the calibration and 
waveform-modeling errors. If these errors can be con- 
sidered uncorrelated, then r) c = 1/ V2 is the appropri- 
ate choice. Waveform- modeling errors, 8h m , are uncor- 
related with the calibration errors, Shu, if (Sh m \Sh R ) 2 <C 
(6h m \5h m }{Sh R \5h R ). To justify taking r\ c = l/y/2, this 
uncorrelated condition should be satisfied for all times 
when data is being collected, and for all the different er- 
rors associated with the various model waveforms that 
will be used in the subsequent data analysis. While cali- 
bration and modeling errors come from completely unre- 
lated sources, it may be difficult to completely rule out 
the possibility of significant correlations. For example, it 
is easy to imagine both calibration and modeling errors 
that are proportional to the exact waveform h e . When 
correlations can not be ruled out, the somewhat stricter 
condition, r\ c — i is required. The calibration accuracy 
needed to ensure no loss of scientific information for the 
strongest sources is generally much greater than the ac- 
curacy needed for detection. The weak dependence of the 
detection standard on rj c has therefore been ignored here, 
and the ideal-detector version of the standards, Eq. (j4}, 
are used in this paper. 



Appendix B: Compound Accuracy Standards 

A somewhat more general class of time-domain accu- 
racy standards are derived in this Appendix using the 
methods developed in SecHU These additional standards 
are somewhat less restrictive for detection than Eq. (|TT)1) , 
but they are also more complex. The idea is to con- 
struct a new measure of the waveform error that uses a 
combination of the time-domain norms used in Eq. (|16p . 
This is done by splitting the integral that appears in the 
basic inequality, Eq. (|12[) . into parts proportional to (i 
and 1 — (j, respectively (where fi is an arbitrary splitting 
factor that satisfies < fi < 1), and then applying the 
constant-noise approximation using different values of k 
for each part. The result is a compound version of the 
basic inequality: 



(6h m \6h r , 



< 



^J oc (2nfr k \dh m (f)\ 2 df 

min[(27r/)2*S n (/)] 
4(1 - M ) f™(2Trf) 2k '\5h m (f)\ 2 df 



rmn[(27r/)^' <?„(/)] ' (B1) 

The right side of this inequality can be re- written in terms 
of the time-domain error measures £ k of Eq. (|17[) . and the 
quantities of Eq. (TTg]) : 



(Sh m \Sh m ) < 



HE* , (1- M )4 



+ 



CI 



(B2) 



This basic inequality can be transformed into a more 
compact and useful form by replacing the splitting pa- 
rameter jj, with a new parameter er, defined by 



/' 



2 ■ 



(B3) 



The parameter /i satisfies the required inequality, < 
fi < 1, whenever < cr < 1. The basic inequality then 
becomes, 



(Sh m \5h m ) < 



c%M' 



(B4) 



where £k,k' (f) is a composite measure of the waveform 
error, defined by 

£l k ,(a)=a£ 2 k + (l-a)£ 2 k ,, (B5) 

and Ck,k'{p) is the composite analog of Ck, defined by 

Cl k ,{a)=aCi + {l-a)Cl. (B6) 



This new form of the basic inequality, Eq. (|B4[) . is true 
for any a. It can be used to convert the fundamental 
frequency-domain standards, Eqs. ([3]) and @, into new 
time-domain waveform accuracy standards: 



£k,k'(&) < Ck,k>{o) 



for measurement, and 



£k,k'(cr) < Ck,k'(<y)v2e 1. 



(B7) 



(B8) 



for detection. The waveform error measure £kk' W) is 
positive definite for all a in the range < a < 1. So 
these new time-domain waveform accuracy standards, 
Eqs. (|B7| and (|B8|) . are sufficient to enforce the funda- 
mental accuracy standards for any a in this range. In 
the extreme cases a = and a = 1, £k,k'{&) reduces to 
one of the single k waveform error measures and these 
standards reduce to special cases of Eqs. (|15|) and (|16|) . 
For intermediate values < a < 1 however, these time- 
domain waveform accuracy standards are completely new 
and distinct. 

Figure 0] illustrates Co^(a) = y / ctCq + (1 — o)C\ as a 
function of mass for the waveform and noise curve used in 
the previous examples. These graphs show that Co,i(c) 
reduces to C Oi i(0) = C\ or C ,i(l) = C in the limit- 
ing cases a = or a = 1, respectively. These graphs 
also show that none of the Co,i(cr) is larger than one 
of the limiting curves, a = or a = 1, for any value of 
A/. Consequently the new compound accuracy standards 
for measurement, Eq. (|B7j) . are generally more restric- 
tive than the simpler single k standard for measurement, 
Eq. (|15p . Therefore, the new compound accuracy stan- 
dards are probably only useful for detection. 

For detection, Eq. (|B8|) determines the maximum al- 
lowed sizes of the error measures £k,k'(&)- Since these 
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M/M, 




FIG. 4: Curves show C ,i(cr) = ^/o-Cg(M) + (1 - a)C7?(M) 
as a function of mass for several values of the parameter a. 



'0 



FIG. 5: Illustrates the allowed regions of the two-dimensional 
space of model waveform errors, (£o,£i), satisfying various 
versions of the accuracy standards discussed in the text. 



measures are independent of the mass of the waveform, 
it follows that these detection standards will be satis- 
fied for the full range of masses only if they are satisfied 
for Cfc.fc'(cr) = minjif Cfe,fe'(cr, M), where the minimiza- 
tion is over the relevant range of masses. Figure [4] shows 
that these minima depend on the parameter a. Perhaps 
the simplest of these compound error measures is the one 
with a = 0.5, in which the two single k error measures, £ & 
and £y , are weighted equally. For the waveforms used in 
this example, minM Co, i (0.5, M) « 0.314, so the detec- 
tion standard becomes £o,i(0.5) < 0.314y'2e m ax. Since 
the factor 0.314 is a little smaller than the factor 0.331 
that appears in Eq. (|19[) . this particular standard is con- 
sidered a little more restrictive. A deeper analysis below 
however shows that this compound standard is more re- 
strictive in some cases but less restrictive in others com- 
pared to Eq. (fig|) . 

There does exist, however, an unconditionally less re- 
strictive compound error standard for a particular op- 
timal choice of the parameter a = <7 op t. Figure 2] 
shows that all the different Co,i(cr) curves intersect at 
the single point M = M int where Co (Mint) = Ci(M; nt ). 
Thus the largest of the minM Co i i(cr, M) can never be 
larger than Co(Afj nt ). The optimal value cr pt, for which 
minM Co,i(<7o P t, M) = Co (Mi n t), can be determined by 
enforcing the extremum condition, dCg i(<r, M)/dM = 0, 
at the point where M = Mint- The result is 

_dC\(dCl dC*Y\ 
aopt -dM{dM-dM) ~°' 284 (B9) 

for the waveform and noise curve example discussed in 
Sec. HI The optimal version of the new compound time- 



domain accuracy standard for detection then becomes: 

£ ,i(0.284)< 0.331^2^. (B10) 

This new version of the detection standard uses the value 
Co.i ~ 0.331, so in this sense it is no more or less restric- 
tive than the previous standard, Eq. (|19|). However as 
the discussion below shows, a larger portion of the two- 
dimensional space of error measures, (£q,£i), is allowed 
using this new standard. So in this sense it is strictly less 
restrictive than the detection accuracy standard given in 
Eq. CE3)- 

To help understand the relationship between these var- 
ious new standards, Fig. [5] illustrates the regions in the 
two-dimensional space of error measures, (£o,£±), al- 
lowed by each version. The original time-domain stan- 
dard, Eq. (|10[) with Co sa 0.016, restricts £q to the region 
left of the vertical dashed line in this figure. The simplest 
version of the new standard for detection, the k = 1 ver- 
sion of Eq. (fT6]) with C\ « 0.239, restricts £\ to the region 
below the dotted horizontal line. A slightly more compli- 
cated standard, Eq. (flU|) with Co = C\ = 0.331, restricts 
both £q and £\ to the region inside the solid-line square. 
The more complicated compound time-domain standards 
presented in this Appendix restrict both £q and £\ to the 
region inside an ellipse defined in Eq. (|B8[) . The simplest 
of these uses the compound error measure £0,1 (cr) defined 
in Eq. (|B5|) with a = 0.5. This limits points to the region 
inside the dash-dotted circle defined by Co,i (0.5) w 0.314. 
The least restrictive standard of this type uses the error 
measure £ ,i (0.284) with C ,i (0.284) 0.331. It limits 
points to the region inside the dash-double-dotted ellipse 
shown in this figure. 

Each of these versions of the waveform accuracy stan- 
dard for detection is sufficient to guarantee the original 
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frequency-domain standards arc satisfied. So any point, 
(£q,£i), in the two-dimensional space of waveform er- 
rors, allowed by any of these conditions is an acceptable 
standard satisfying point. The most general accuracy 
standard of this type could be obtained by taking the 
union of the allowed regions inside each of the ellipses 
defined by Eq. flB8|) with C ,i(a) = min M C ,i(ct, M). 



The boundary of this region qualitatively has the shape 
of a hyperbola that asymptotes to the vertical dashed 
line and the horizontal dotted line in Fig. [5] This re- 
gion is rather complicated to specify exactly, however, so 
this most general version of these compound standards 
is probably impractical for widespread use. 
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